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Abstract 

The transverse momentum and rapidity distributions of net protons and negatively 
charged hadrons have been measured for minimum bias proton-nucleus and deuteron- 
gold interactions, as well as central oxygen-gold and sulphur-nucleus collisions at 200 
GeV per nucleon. The rapidity density of net protons at midrapidity in central nucleus- 
nucleus collisions increases both with target mass for sulphur projectiles and with the 
projectile mass for a gold target. The shape of the rapidity distributions of net pro- 
tons forward of midrapidity for d+Au and central S+Au collisions is similar. The average 
rapidity loss is larger than 2 units of rapidity for reactions with the gold target. The trans- 
verse momentum spectra of net protons for all reactions can be described by a thermal 
distribution with 'temperatures' between 145 ±11 MeV (p+S interactions) and 244 ± 43 
MeV (central S+Au collisions). The multiplicity of negatively charged hadrons increases 
with the mass of the colliding system. The shape of the transverse momentum spectra 
of negatively charged hadrons changes from minimum bias p-|-p and p+S interactions 
to p+Au and central nucleus-nucleus collisions. The mean transverse momentum is al- 
most constant in the vicinity of midrapidity and shows little variation with the target 
and projectile masses. The average number of produced negatively charged hadrons per 
participant baryon increases slightly from p+p, p+A to central S+S,Ag collisions. 
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1. Introduction 

In order to investigate and understand the dynamics of relativistic nucleus-nucleus col- 
lisions, it is important to have information on the rapidity and transverse momentum 
distributions of nucleons participating in the interaction and of particles produced in 
the collisions, all as a function of the size of the system and of the impact parameter. 
The differences of particle distributions from proton-proton, proton-nucleus and finally 
nucleus-nucleus interactions may provide the key to understanding any non-hadronic ef- 
fects that may appear in central nucleus-nucleus collisions. Ultimately these distributions 
might be used to determine the behavior of strongly interacting matter at high densities 
and temperatures. 

New results from the measurement of collisions of very heavy nuclei at the BNL AGS 
and CERN SPS [ ^ § can best be understood by comparison to simpler systems of lighter 
projectiles and targets. We report results from the NA35 experiment at the CERN SPS 
on the rapidity and transverse momentum distributions of negatively charged hadrons 
h~ and participant protons, measured in proton-nucleus, deuteron-gold, oxygen-gold 
and sulphur-nucleus collisions at 200 GeV per nucleon. A small subset of the present 
data (S+S collisions) has already been published [ |^ and has been compared to various 
theoretical predictions ([^, |^, Q). The transverse momentum and rapidity distributions 
provide information on the degree of stopping, thermalization, expansion and flow effects 
in the collisions. Furthermore, the participant proton rapidity distributions may be used 
to determine the baryon density or the baryo-chemical potential as a function of rapidity, 
which constrains model calculations of the dynamics of particle production. The present 
study is complementary to transverse and forward energy measurements already reported 
from NA35 [ which provided global information on the collision dynamics. It is also 
supplementary to previous strange particle production measurements [ ||, ^ |10|, |Tl 



which are the basis for the determination of the chemical freeze-out conditions. All of 
this information helps to obtain a consistent picture of the entire reaction process. In 
particular, it will facilitate detailed comparisons with model calculations (e.g. [ ^ |T^, 

m B 0)- 
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2. Experiment NA35 

Experiment NA35 at tlie CERN SPS studies collisions of p, d, ^^O and '^^S projectiles of 
200 GeV per nucleon incident energy (yj/viv=19.4 GeV) with nuclear targets. The main 
detectors are two large- volume tracking devices: a streamer chamber (SC) inside a 1.5 
Tesla vertex magnet (VTM) and a time projection chamber (TPC) positioned downstream 
of the magnet. Both detectors record the space trajectories of charged particles from 
which the particle momenta are derived. The acceptances of the detectors for negatively 
charged hadrons are complementary to each other and cover the full phase space with 
some overlap. Central events are usually selected by the absence of projectile spectators 
in the veto calorimeter placed in the beam line downstream of all detectors (see e.g. [ 

S 0, 0)- 

The objective of the experiment is to measure the hadronic final state in nucleus- 
nucleus collisions. In addition to negatively and positively charged particles, weak decays 
of charged and neutral particles can be detected in the streamer chamber. The analysis 
of the streamer chamber pictures has been fully automated [ 0, which not only resulted 
in an increase of the measuring rate by an order of magnitude but also allowed for a com- 
plete simulation of the detector based on the GEANT-package. The results discussed in 
the following are mainly based on the automated analysis but traditional measurements 
have been included where appropriate. The TPC identifies charged particles by multiple 
sampling of the specific ionisation with a dE/dx resolution of about 5-6%. It covers the 
momentum range of the relativistic rise of the energy loss. This makes particle identifica- 
tion using statistical methods possible. Combining both detectors, spectra of negatively 
charged hadrons i^^, , K~ , net protons {p — p), A, A and p are obtained. Thus 
information on practically all baryons and mesons in the final hadronic state can be ac- 
cessed with large phase space coverage. This paper will focus on rapidity and transverse 
momentum spectra of net (mainly participating) protons and negatively charged hadrons 
(primarily pions). 

In order to accomplish large phase space coverage in collisions with very different par- 
ticle multiplicities (from p+S to S+Au) and track densities in the laboratory, data were 
taken in three different experimental setups. Proton- and the oxygen-induced reactions 
were studied using the (standard) streamer chamber configuration (SETUP I) described 
in [ ^, 16, 17]. The nuclear target was placed 11 cm upstream of the 2.0 x 1.2 x 0.72 
m^ streamer chamber (external target) which is the main tracking device inside the 1.5 T 
magnetic field of the superconducting vertex magnet. In addition p+S interactions were 
recorded with a target placed inside the streamer chamber (internal target). Charged 
particles created tracks of streamers of about 1 mm diameter and 1 cm length along 
their path through the sensitive volume. Three cameras, each equipped with a two-stage 
magnetically focusing image intensifier of about 2500-fold gain, recorded the events with 
a demagnification of 30 on 70 mm film. For the ^^S induced reactions the track den- 
sity in the central region along the beam direction {Qiab < 5°) is so high that individual 
tracks cannot be separated due to the projective readout of the streamer chamber by 
stereoscopic cameras. Therefore, the streamer chamber geometry was modified by the 
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insertion of an opaque foil of a trapezoidal shape close to the beam axis reducing the 
projected track density in that region. In addition, another opaque foil screened the lower 
half of the streamer chamber in order to further reduce the track density. This modified 
configuration is called SETUP II. 




Figure 1. Experimental setup of experiment NA35 at the CERN SPS. The target is placed 
either in front of the streamer chamber (SETUP I and II) or in front of the sweeper magnet 
SM (SETUP III). 



In order to measure forward-going strange particles and participant protons in the 
projectile hemisphere for the heaviest system (S+Au), the target and a sweeper magnet 
(SM) were installed upstream of the standard SETUP I configuration as shown in Fig. ^ 
(called SETUP III). The sweeper magnet has a vertical entrance opening of ±2.5° for 
particles emerging from the target placed 6.3 cm upstream of this magnet. Due to its 
bending power of about 4.5 Tm all low momentum charged particles are swept out of 
the acceptance of the second magnet thus reducing the track density in the streamer 
chamber. This configuration not only enables measurement of strange particle decays in 
the projectile hemisphere but also covers participant protons beyond midrapidity (?/cm = 
3.0) up to beam rapidity {i/proj = 6.0). This can be seen in Fig. ^ where the population of 
protons in transverse momentum pt and rapidity y is presented. When sulphur projectiles 
are used the proton acceptance can be extended up to beam rapidity only in SETUP III. 

A second large-volume tracking device, a time projection chamber with a sensitive 
volume of 240 x 125 x 112 cm^ [ Il8[ located downstream behind the analysing magnet 



was installed in 1990 and utilized in the ^^S-beam running for all three setups. It was also 
used in the deuteron-nucleus studies in the SETUP I configuration. The TPC detects 
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charged particles by measuring the ionisation of a particle, traversing the TPC, up to 60 
times. For a primary track emerging from the production vertex the momentum of the 
particle can be calculated from the direction and position of the track inside the TPC. 
The TPC was positioned for most of the time such that it detected magnetically stiff 
negatively charged particles above midrapidity. The acceptance for negatively charged 
hadrons (pions) is illustrated in Fig. |^. Both tracking detectors combined cover almost all 
momentum space available for hadron production at SPS energy. In addition, the TPC 
identified charged particles by multiple sampling of the specific ionisation. This capability 
was used to measure antiproton spectra [ |l9l and to distinguish electrons from negatively 
charged hadrons. 
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Figure 2. Acceptance of the streamer chamber in configuration SETUP II (left) and 
SETUP III (right) for protons, illustrated by the phase space population in a typical 
heavy ion collision. Note the forward acceptance of SETUP III. 
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Figure 3. Acceptance of the two large-volume tracking detectors, the streamer cham- 
ber (left) and TPC (right) for SETUP II, illustrated by the phase-space population of 
negatively charged hadrons in a typical heavy ion collision. 
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3. Data sets 

NA35 has studied a variety of collision systems, ranging from minimum bias proton- 
sulphur interactions to central sulphur-gold reactions. The central-collision fractions of 
the inelastic cross sections of the nucleus-nucleus collisions are typically the upper few 
percent of the inelastic cross section for the heavy ion induced reactions (see Table |l]). 
Central events are selected either by a minimal energy deposition in the projectile frag- 
mentation domain at forward angles [Qiab < 0.3°) which is covered by the veto calorimeter 
(FET trigger), or by a large transverse energy deposition in the midrapidity calorimeters 
(TET trigger). In the case of p-|-S interactions a trigger was derived from a scintilla- 
tion counter placed downstream of the streamer chamber which covered the projectile 
fragmentation domain. Its veto threshold was set just below the signal of the projectile 
particle and it selected essentially all inelastic interactions. An additional off-line selec- 
tion of events with a charged particle multiplicity larger than five was applied in order to 
exclude contaminations due to non-target interactions in low multiplicity events. 

Table 1 

Characteristics of the data sets: colliding system, experimental configuration, number 
of events, trigger type, the selected fraction of the inelastic cross section, acceptance 
in rapidity for negatively charged hadrons and for protons. (The trigger for the d+Au 
interactions marked by (*) selects collisions where the deuteron projectile fully overlaps 
with the gold target.) 



Reaction 


Setup 


Number 
of events 


Trigger 




Acceptance 
h- 


Acceptance 
protons 


P+S 


SETUP I 




scint. veto 

h+- > 5 










internal target 


735 




63.5 


-2<y^<4 


-1 < 1/p < 3 




external target 


339 




63.5 


0.6 <y^<7 


0.2 <yp<6 


p-t-Au 


SETUP I 


1603 


FET 


67 


0.6 <y^<6 


0.2 <yp<6 


d+Au 


SETUP I 


2457 


FET* 


43 


0.4 <y^<6 


0.2 <yp<6 




TPC 


9484 


FET* 


43 




2.8 <yp< AA 


0+Au 


SETUP I 


565 


FET 


1.8 


0.6 <y^< 3.8 


0.2 < ?/p < 3 


s+s 


SETUP II 


2650 


FET 


2.9 


0.6 <y^ < 3.8 


0.2 < < 3 




TPC 


23500 


FET 


3.5 


3.2 <y^ < 5.4 




S+Ag 


SETUP II 
TPC 


13730 
18000 


FET 
FET 


3.2 
3.2 


0.6 <y^< 3.8 
3.2 <y^< 5.4 


0.2 < ?/p < 3 


S+Au 


SETUP I 


331 


TET 


1.3 


0.8 <yn< 3.6 






SETUP III 


3142 


FET 


6.3 


3.8 <y^ < 6.0 


2.8 < yp < 6 




TPC 


25700 


FET 


6.3 


3.2 <y^< 5.4 
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4. Data Analysis 



The operation of tlie streamer ciiamber (SC) is described elsewhere [ |I6|, 0, the data 
analysis will briefly be described in the following. Streamer chamber pictures were scanned 
and measured using two different approaches: using the traditional fllm measurement by 
human operators or applying image processing and pattern recognition techniques to 
digitized streamer chamber data. Most of the streamer chamber data were analyzed 
automatically by the Frankfurt Image Processing System which is described in detail 
elsewhere [ |16[ . The three perspective views of a streamer chamber event were recorded on 



70 mm fllm and later digitized by a high resolution CCD-linescan camera. These digitized 
images were submitted to a series of image processing and pattern recognition programs. 
Tracks matched in the three stereoscopic views were then reconstructed in space. The 
particle momenta and the production vertex were obtained by fltting a trajectory to 
the measured points. The low multiplicity data obtained in SETUP I were analyzed on 
traditional projection tables by human operators. 

Since the time projection chamber (TPC) was new to NA35 and relatively new for 



heavy ion applications [ |T8|, we will describe its operation briefly. A charged particle 
which traverses the sensitive volume of the TPC, ionizes gas molecules and produces 
electrons along its trajectory. Guided by a homogeneous electric drift fleld, the electrons 
drift towards the read-out MWPCs of the TPC. A fleld gradient at the sense wires leads 
to a multiplication of the drifted electrons by 5 ■ 10^. The resulting space charge induces 
a signal onto the segmented cathode plane (pads). The centroid of these induced signals 
(clusters) can be determined with an accuracy of o" = 500 fim in pad direction and of 
cr = 310 /im in drift (time) direction. Pattern recognition is done by a simple track 



follower, which associates clusters with a straight line [ |20| . The particle momenta are 



determined by assuming that the particles emerge from a known production vertex. 

4.1. Track recognition and reconstruction efficiency 

The momentum resolution achieved with the streamer chamber as a continuous tracking 
device is approximately 3 MeV/c {dp/p = 0.15%) for a typical pion with momentum of 
2 GeV/c. The momenta in the TPC are typically higher than in the streamer chamber, 
since the TPC is located at more forward angles in the laboratory. The TPC resolution 
for 10 GeV/c tracks is approximately 1%. 

In the streamer chamber track reconstruction losses due to tracks which overlap in the 
individual camera projections as well as to two-track resolution, pattern recognition and 
space reconstruction inefficiencies can be examined by either of two methods. The first 
method relies on human visual scanning of film data and multiple independent scans of the 
same data. The second method is a Monte Carlo simulation of digitally processed images. 
By this method heavy ion events are simulated using the FRITIOF event generator code [ 
pT|. The track trajectories are modeled through the NA35 experimental set-up using the 



GEANT detector simulation package [|2^. For each track the streamers are calculated and 
then projected through the imaging system onto film. The film response and digitisation 
process are then simulated. The resulting digital event pictures are submitted to the 
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entire pattern recognition and track reconstruction process. This procedure provides 
the efficiency as a function of the two kinematical variables rapidity y and transverse 
momentum p^- Both methods give an overall efficiency of about 95-97% for the highest 
multiplicities in our acceptance. 

For the TPC the reconstruction efficiency was determined by a Monte-Carlo simulation 
using the VENUS event generator [ |12|], the GEANT-program and a simulation of the 
TPC response [ The efficiency was found to be 90-95%, approximately independent 
of y and px within the geometrical acceptance. For the particle spectra the inefficiencies 
were corrected function of y and px- 

A comparison of rapidity densities of negatively charged hadrons obtained by the 
streamer chamber analysis of SETUP II and SETUP III data and by the TPC analy- 
sis is shown in Fig. ^. For the symmetric system S+S streamer chamber data below 
midrapidity can be reflected about ycm = 3.0 and then compared to TPC data above 
midrapidity as shown in Fig. |^. The TPC data below y = 4.0 and above y = 4.6 have 
been extrapolated to low and high pt respectively by an exponential in l/pr ■ dn/dpT- 
The comparisons in both figures demonstrate a satisfactory agreement (within several %). 
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Figure 4. Rapidity distributions of negatively charged hadrons (pion mass assumed) 
produced in central S+Au collisions, measured in the streamer chamber (configurations 
SETUP I and SETUP III) and in the TPC. 
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Figure 5. Rapidity distributions of negatively charged hadrons produced in central S+S 
collisions, measured in the streamer chamber (SETUP II) and reflected around i/cm = 3, 
and in the TPC. 



4.2. Particle identification 

We define negatively charged hadrons h~ as the sum of vr", K' and p. Their aver- 
age multiplicities in central S+S collisions are estimated to be about 90:7:1 [ |TT], |T9[. 
Corrections for electron contamination were made on a statistical basis (see next chap- 
ter). The rapidity of negatively charged hadrons was calculated assuming the pion mass, 
throughout; this implies that the K~ and p contributions to h~ rapidity distributions 
appear with an upward shift relative to their true position. A slight asymmetric widening 
of the rapidity distribution results. Since the streamer chamber has very little particle 
identification capability for stable particles, the net proton distribution was deduced from 
a measurement of the charge excess 

where the rapidity is calculated assuming the proton mass. This is a good approximation 
for all rapidities in collisions of isoscalar nuclei (equal number of neutrons and protons) 
assuming local isospin conservation. Even for asymmetric systems deviations are expected 
to be small. Equating net protons to the charge-excess is subject to corrections which 
will be discussed below. 
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4.3. Corrections 

All data in this paper are corrected for geometric acceptance and reconstruction inef- 
ficiencies, for and e~ misidentified as hadrons, for hadrons resulting from the weak 
decays of Kg, hyperons and for hadrons produced in secondary hadron- nucleus interac- 
tions. 
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Figure 6. Transverse momentum spectrum of conversion electrons in central S+S collisions 
measured and identified in the TPC in the pion rapidity range 4.0 < y < 4.4 compared 
to a calculated spectrum based on the measurement of negatively charged hadrons (see 
text for details). The vertical scale is given in (GeV/c)~^. 



One source of contamination in the charged hadron spectra is the 27-decay of the 
abundantly produced 7r° inside the target and subsequent conversion of 7s into e'*"e^-pairs. 
The spectral distributions of contaminant electrons were obtained from a Monte Carlo 
calculation using the experimental negative particle distributions and assuming equal 7r° 
and 7r~ meson multiplicities and spectra. An alternative method is to use the vr" spectra 
obtained from an event generator (FRITIOF), with parameters tuned to reproduce the 
measured spectra, and a GEANT simulation for the subsequent decays and conversions. 
The latter method was employed at high rapidities, where the data based Monte Carlo 
may suffer from reconstruction inefficiencies at high momenta. The methods give similar 
results, which are confirmed by identified electron spectra that are measured in the TPC 
and displayed in Fig. |[ The electron contribution to the charged particle spectra is 
subtracted for each (?/,Pt) bin. 

The contamination of the negatively charged hadron spectra by pions produced in 
secondary interactions primarily in the target, and in weak decays is corrected for on the 
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basis of GEANT simulations (see e.g. [|2^). The parameters of the event generators were 
tuned to reproduce the measured strange particle yields and spectra. 

The contributions of the various contaminations as functions of pt and y for a central 
S+Ag collision are shown in Figs. and ^ respectively. The contaminations are largest 
at small pt- The electrons dominate below transverse momenta of 150 MeV/c. The first 
Pt bin (50 MeV/c) especially is strongly affected by the electron correction, the amount 
of the contamination can be up to 50%. Above 150 MeV/c the weak decay products 
represent the main contamination to the spectrum. For larger transverse momenta the 
various contributions get smaller. The rapidity dependence of the decay product contami- 
nation follows the distribution of the primaries, while electrons and secondary interactions 
contribute more strongly at low rapidities (Fig. H). The overall systematic error on the 
spectra is estimated to be less than 5% for px above 50 MeV/c. 

All data - spectra of both negatively and positively charged hadrons - were corrected 
as described above. In addition, when deriving the net proton spectra from the charge 
excess, a correction is applied for the excess of the relative to the K~ yield. Since these 
corrections are based on measurements of the same experiment over a wide acceptance, 
the corrections are model independent. Nevertheless, these corrections become large at 
low transverse momenta, particularly at low rapidities, and therefore a low-pr cut was 
applied to the net proton spectra. This cut also takes care of possible differences in the 
TT"'" and 7r~ yields for non-isoscalar systems of colliding nuclei. In order to obtain rapidity 
densities the transverse momentum spectrum was extrapolated into the part below the 
cut using a thermal fit to the data. This extrapolation contributes approximately 20% 
(assuming a temperature of about 200 MeV and a px-cut of 300 MeV/c) to the cross 
section. 

The validity of equating participant protons to the charge-excess even in asymmetric 
systems is demonstrated in Fig. ^ Particles taken from a simulation (FRITIOF) were 
passed to a GEANT-program of the experimental setup. After all corrections described 
above the net proton distribution which was deduced from the charge excess agrees very 
well with the input distribution of net protons within the acceptance 0.2 < y < 3.0. The 
overall systematic error on the participant proton spectra is estimated to be less than 12%. 
Note that this conclusion does not rest on the implication that the FRITIOF prediction 
in Fig. 1^ gives a good account of the data (see sect. 6.1.3). 
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Figure 7. Composition of a measured transverse momentum spectrum of negatively 
charged particles, based on a tuned event generator (FRITIOF) and a GEANT simu- 
lation. The vertical scale is given in (GeV/c)^^. 




Figure 8. Composition of a measured rapidity spectrum of negatively charged particles, 
based on a tuned event generator (FRITIOF) and a GEANT simulation. 




Figure 9. Simulated rapidity distributions of net protons {p — p) and of reconstructed net 
protons deduced from the charge excess after all corrections for central S+Ag coUisions. 
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5. Net Baryons 

In a simple picture of high energy nucleus-nucleus collisions one distinguishes two 
classes of nucleons: spectators and participants. Spectator nucleons are defined as those 
which did not suffer collisions and therefore have their momenta located inside the Fermi 
spheres of the interacting nuclei. Participant nucleons are those which suffer collisions 
that affect their momentum sufficiently to remove them from the Fermi sphere of the 
incoming nuclei. Some of them may be converted into strange and multistrange baryons. 
Integrating the rapidity distributions of net baryons (baryons minus antibaryons) between 
y = 0.2 and y = 5.8 yields an estimate of the total number of participant baryons. This 
estimate is only accurate in the case of central collisions of identical nuclei (S+S), where 
the number of spectators is small compared to the number of participants so that the 
result of the technical procedure to separate spectators and participants depends only 
weakly on the cuts used. For asymmetric systems the net baryon number as calculated 
in this paper should be treated as a crude approximation of the number of participants 
since cascading processes in the target spectator matter may invalidate the definition of 
a spectator. 

5.1. Rapidity distributions 

Rapidity distributions of net protons were obtained from the charge excess with the 
associated corrections discussed above. Since these corrections become large at low trans- 
verse momenta, a cut was applied at 300 MeV/c (for S+A collisions at 400 MeV/c). The 
missing low momentum part was extrapolated on the basis of a thermal fit in order to 
obtain rapidity densities (see chapter 5.2). In comparing rapidity distributions for the 
various systems one has to keep in mind that different forward energy triggers may affect 
the spectra differently. 

5.1.1. Proton-Nucleus Collisions 

The rapidity distributions of net protons for minimum bias p+S and p+Au interactions 
are shown in Fig. 0. While the rapidity densities around midrapidity {ycm = 3) for both 
reactions are not very different, clear differences are observed for rapidities below 1.2 and 
larger than 4.4. In p+Au collisions more target nucleons participate in the reactions and 
are therefore shifted up to one unit in rapidity farther than in p+S collisions. For the 
projectile the gold nucleus looks black, i.e. the probabihty of the projectile to traverse 
the target nucleus without losing any or little energy is small. For p+S interactions on 
the other hand, one can perceive a clear projectile proton peak at rapidity ~ 4.7. 
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Figure 10. Rapidity distributions of net protons {p—p) for minimum bias p+S and p+Au 
collisions at 200 GeV/nucleon. 



5.1.2. Nucleus-Gold Collisions 

In order to study the dependence of nuclear stopping on the mass of the projectile in 
the case of incidence on a heavy target nucleus, net proton rapidity distributions were 
measured for various 'full overlap' collisions: central d+Au, central 0+Au and central 
S+Au reactions. These data are displayed in Fig. |Tl|. 0+Au collisions were measured in 
the streamer chamber for rapidities below midrapidity, y < i/cm, SETUP 1. The S+Au data 
were measured in the streamer chamber configuration SETUP 111 at rapidities y > ycm- 
The rapidity density at midrapidity, y = ycm, is observed to increase with projectile mass. 
A comparison of the rapidity distribution for S+Au collisions with the distributions for 
d+Au and 0+Au, multiplied by 16 and 2 (the inverse ratios of the projectile nucleons 
relative to the S+Au reaction), respectively, is displayed in Fig. |12|. For the 0+Au data 
at a beam energy of 60 GeV/nucleon the rapidity axis is scaled by the ratio of the beam 
rapidities {yproj (200 GeV/nucleon) / yproj (60 GeV/nucleon) = 6.0/4.8) prior to plotting 
in Fig. |l^ [ The resultant rapidity distributions are similar for y > ycm- This 

implies that the stopping of projectile nucleons is rather independent of the size of the 
projectile nucleus and its energy in collisions with a heavy target nucleus. 
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Figure 11. Rapidity distributions of net protons (p — p) for central d+Au, central 0+Au 
and central S+Au collisions at 200 GeV/nucleon. 
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Figure 12. Scaled rapidity distributions of net protons for central 16- (d+Au), central 
2-(0+Au) and central S+Au collisions. For the case of the 60 GeV/nuclcon data the 
rapidities are scaled by the ratio of the beam rapidity relative to that of 200 GeV/nucleon. 
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5.1.3. Sulphur-Nucleus Collisions 

To investigate tlie effect of clianging tlie nuclear tliickness of tlie target, rapidity distri- 
butions of net protons were measured for central S+S, S+Ag and S+Au collisions. The 
data are shown in Fig. The S+S and S+Ag data were measured in the SETUP 11 for 
y < Ucm and S+Au data in the SETUP 111 for y > i/cm- For the symmetric system S+S, 
the data points reflected at i/cm are also shown. The shape of the net proton rapidity 
distribution changes from being relatively flat for the light symmetric system (S+S) to 
a monotonic decrease with rapidity for the asymmetric systems. The rapidity density at 
y = Vcm increases for heavier projectiles, and is a factor of about two larger for S+Au 
relative to S+S collisions. There is an excess of net protons in the forward direction for 
the lighter S+S reaction relative to the S+Au reaction. Thus, it is evident that the S 
target nucleus is less effective than the Au nucleus in stopping (slowing down) incident 
projectile nucleons. For the S+Au system, projectile protons are shifted further backward 
in rapidity than for the S+S system. The rapidity densities near midrapidity agree with 
the values obtained by the NA44 experiment for identified protons and antiprotons as 
shown in Table 0. 
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Figure 13. Rapidity distributions of net protons {p — p) for various central S-nucleus 
collisions at 200 GeV/nucleon. 



5.1.4. Rapidity Shift 

Before the collision the nucleon rapidity distribution is concentrated in peaks at the tar- 
get rapidity y = and at beam rapidity y = yproj, with the peak width given by the Fermi 
momentum inside the nuclei {Ay ^ 0.2). After the collision, target and projectile nucle- 
ons are shifted towards midrapidity, with some nucleons converted into (multi) strange 



20 



Table 2 

Rapidity densities of net protons {p 



p) near midrapidity {y ~ 2.8) in central S+S and 



S+Au(Pb) collisions [ |7 





dn/dy {p — p) 


Reaction 


this paper 


NA44 [ |23] 


S+S 


3.2 ± 0.4 


3.2 ± 1.5 


S+Au(Pb) 


12 ± 2 


12.4 ± 2.1 



baryons. The latter must also be counted in the evaluation of rapidity loss. The rapidity 
distributions of net hyperons (A — A) for minimum bias p+S and p+Au collisions are 
displayed in Fig. [1^. Those for central collisions of S+S, S+Ag, and S+Au are displayed 
in Fig. |15[ The trends in the rapidity distributions for net hyperons as a function of target 
nucleus are similar to those for net protons described in the previous section. 

In order to quantify the amount of stopping in symmetric and asymmetric collisions, 
the average rapidity loss (or mean rapidity shift) {5y) = yproj — {y) is defined, where yproj 
is the incoming projectile rapidity and {y) is the average net baryon rapidity after the 
collision. For symmetric systems the average rapidity was calculated by averaging the 
rapidity of the net protons in the interval from ycm = 3 to yproj — 0.2. For the asymmetric 
systems {y) is defined as the average rapidity of the net proton dn/dy distribution in 
the projectile hemisphere from yproj — 0.2 to a rapidity such that the integral under the 
distribution equals the number of net projectile protons If the target and projectile 
nucleons mix, the rapidity loss defined in this manner is a lower limit. 

The average rapidity losses {5y) for collisions of various projectiles with sulphur and 
gold targets are displayed in Fig. |l^. For each target {5y) is approximately independent 
of projectile mass (but may also be construed as being consistent with a slight increase 
with projectile mass). The most obvious trend is a shift from {5y) = 1.2 - 1.6 for sulphur 
target interactions to {6y) = 2.0 - 2.3 for the gold target interactions. 
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Figure 14. Rapidity distributions of net hyperons (A — A) for minimum bias p+S (1.0 < 
y < 5.0) and p+Au (1.4 < y < 4.4) collisions at 200 GeV/nucleon. 
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Figure 15. Rapidity distributions of net hyperons (A — A) for central S+S (0.5 < y < 3.0), 
S+Ag (0.5 < y < 3.0) and S+Au (3.0 < y < 5.0) collisions at 200 GeV/nucleon. The 
solid stars are measured S+S data, the open stars are reflected at y = ycm- 
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Figure 16. Average rapidity loss (mean rapidity shift) of net projectile protons for various 
projectiles incident on sulphur and gold targets at 200 GeV per nucleon. 



5.1.5. Number of Net Baryons 

By integrating the rapidity distributions the number of net protons inside the accep- 
tance was obtained and is shown in Table ^. The total number of net protons is defined as 
the net proton number in the final state outside the 'Fermi spheres' of the projectile and 
target nuclei - approximated by selecting a rapidity range 0.2 to 5.8. For systems where 
the acceptance does not cover the rapidity range from 0.2 to 5.8, we extrapolated the 
rapidity distribution. For the O+Au data the extrapolation was based on the observed 
scaling behaviour of 60 GeV/nucleon 0+Au data (see section 5.1.2). The number of net 
protons in S+Ag coUisions at 3.0 < y < 5.8 was estimated to be the average of S+S and 
S+Au data. 

The charge-excess method, corrected for the K^/K~ asymmetry and strange hyperon 
decays, leads to the number of observednet protons. In order to estimate the total number 
of net baryons one needs in addition an estimate of the total number of net hyperons. 
Table § shows the measured yields of net lambdas (A— A) measured in the same experiment 
[ 1^, |lOl . We derived the corresponding correction for the total number of net hyperons 
from the (A — A) distributions using Wroblewski's empirical rule [ = 1.6 ■ A^. 

Multiplying the number of net protons by 2 (in the case of p-induced reactions only in 
the target hemisphere) and adding the corrected number of hyperons gives an estimate of 
the total number of net baryons. These estimates are shown in Table ^. This procedure 
is valid for collisions of isoscalar nuclei (equal number of neutrons and protons) and 
only a crude approximation in all other cases. Only in the case of central collisions of 
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isoscalar nuclei the number of net baryons equates the number of participant nucleons. 
Comparing the measured number of net baryons with the number of wounded nucleons 
using the standard geometrical approach [ Q one finds a discrepancy between these two 
numbers in asymmetric systems. This can be explained by target spectators misidentified 
as participants or additional participants due to cascading effects in spectator matter. 



Table 3 

Number of observed net protons (p — p) in p-nucleus, central deuteron- nucleus and central 
nucleus-nucleus collisions. The values in italics are the result of an extrapolation. See 
text for details. 



p — p 


Reaction 


0.2 <y < 3.0 


3.0 < y < 5.8 


0.2<y< 5.8 


p+S 






2.4 ± 0.2 


p+Au 






5.9 ± 0.2 


d+Au 






10.7 ± 0.2 


0+Au 


42 ± 5 




47 ± 5 


S+S 


10.6 ± 0.8 


10.6 ± 0.8 


21.2 ± 1.3 


S+Ag 


32 ± 2 




43 ± 3 


S+Au 




11.4 ± 0.5 





Table 4 

Number of observed net lambdas (A — A) in p-nucleus and central nucleus-nucleus colli- 
sions. The values in italics are the result of an extrapolation. See text for details. 



A- A 


Reaction 


1.0 < y < 5.0 


1.4 < y < 4.4 


0.2 < y < 5.8 


p+S 


0.19 ± 0.02 




0.23 ± 0.03 


p+Au 




0.20 ± 0.02 


0.28 ± 0.03 




0.5 <y< 3.0 


3.0 < y < 5.5 


0.2<y< 5.8 


S+S 


3.6 ± 0.4 


3.6 ± 0.4 


7.2 ±0.6 


S+Ag 


8.5 ± 0.8 


3.9 ± 0.5 


12.6 ± 1.2 






3.0 < y < 5.0 




S+Au 




3.9 ± 0.4 
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Table 5 

Estimated total number of net baryons < B — B > in p-nucleus and central nucleus- 
nucleus collisions compared to the number of wounded nucleons < W > (the values for 
p+S and p+Au marked by (*) are averaged over all impact parameters). 



Reaction 


<B-B> 
0.2 < y < 5.8 


<W> 
calculated 


p+S 


4.2 ± 0.8 


3.2* 


p+Au 




4.7* 


d+Au 




10 


0+Au 




68 


S+S 


54 ± 3 


53 


S+Ag 


105 ± 12 


92 


S+Au 




113 
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5.2. Transverse momentum distributions of net protons 

The invariant cross sections as a function of the transverse momentum for net protons 
are compatible with the following parametrization of the transverse momentum spectra 
as suggested by thermodynamic models: 

C ■ rriT ■ Ki { 



Pt dpT Ti 



where irtT = yPT + is the transverse mass of particles with rest mass mo, C is a 
normalization factor and Ki is the modified Bessel function of order one. The parameter 
Ti is determined by fitting the experimental data and can be interpreted in thermal models 
without transverse fiow as the thermal freeze-out 'temperature'. 

Figs. p!7|-p!9| show transverse momentum distributions of net protons for various collision 



systems at 200 GeV/nucleon together with fits to Eq. 1. The invariant spectra of protons 
as a function of px for various systems are well described by inverse slope parameters 
('temperatures') of approximately 200 MeV (Table P). The 'temperatures' are observed 
to increase with the mass of the colliding system and are highest at midrapidity. The 
'temperatures' range from T ^ 145 MeV for minimum bias p+S to T ^ 220 — 240 MeV 
for central collisions of S+S, S+Ag and S+Au. Dividing the rapidity range into two bins, 
0.5 < y < 2.0 and 2.0 < y < 3.0, and fitting the spectra separately in these bins results 
in inverse slope parameters that are about 40 MeV higher at midrapidity as compared to 
the values at lower rapidity for all systems. 




Pt [ GeV/c ] 



Figure 17. Transverse momentum distributions of net protons for minimum bias p+S and 
p+Au collisions. The vertical scale is given in (GeV)^^. 
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Figure 18. Transverse momentum distributions of net protons for central d+Au, 0+Au 
and S+Au collisions. The vertical scale is given in (GeV)"'^. 
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Figure 19. Transverse momentum distributions of net protons for central S+S, S+Ag and 
S+Au collisions. The vertical scale is given in (GeV)^^. 
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Fitting a simple exponential 

1 dn mr. , . 

— -^ = C-exp {-—) 2 
Pt dpT J 2 

results in somewhat higher 'temperatures' as shown in Table ^. Also shown are our pub- 
lished results for S+S collisions and results from another experiment. The 'temperature' 
reported by NA35 in [ Q was about 40 MeV lower. Much higher statistics of the recent 
data allowed a more detailed contamination correction as described in Section 4. NA44 
data [ |2^, |3l| measured in a small rapidity interval near midrapidity for less central 
collisions yield somewhat lower inverse slope parameters. 



Table 6 

'Temperatures' obtained from thermal fits (Eq. 1) and fits according to Eq. 2 to the trans- 
verse momentum distributions of participant protons in various collisions. The systematic 
error is about 10 MeV. 





(Eq 


AeV] 
■ 1) 


T2 

(Eq 


AeV] 
•2) 


Reaction 


0.5<y< 3.0 


NA35 


0.5 <y< 3.0 


NA44 


P+S 


145 ± 11 




154 ± 14 




p+Au 


156 ± 4 




163 ± 5 




d+Au 


161 ± 2 




172 ± 5 




0+Au 


204 ± 4 




219 ± 5 




S+S 


226 ± 5 


180 ± 14 [ 1 


235 ± 9 


208 ± 8 [ 


31 


1 


S+Ag 


224 ± 2 




238 ± 2 






3.0<y< 5.0 




3.0 <y< 5.0 




S+Au 


244 ± 43 




276 ± 48 


242 ± 3 [ 
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6. Negatively Charged Hadron Production 

Essentially all of the energy deposited in the reaction volume by the participant nu- 
cleons goes into the production of hadrons. Without extensive particle identification in 
this experiment, only distributions of inclusive charged hadrons can be measured and 
separation of various species (such as tt"*", K"*", p) is not possible. The positive hadron 
distributions are more complicated to interpret than the negative ones, since the contribu- 
tion of and p is a large contamination. Thus, we will concentrate on the distributions 
of negative hadrons which are primarily vr" mesons. The relevant contamination and 
acceptance corrections as described in section 4 were applied. 

6.1. Rapidity Distributions 
6.1.1. Proton-Nucleus Collisions 

The rapidity distributions of negatively charged hadrons produced in p+S and p+Au 
interactions are shown in Fig. For comparison nucleon-nucleon reactions are also 
included. The nucleon-nucleon interactions are defined here as 0.5-[{p+p —>■ TT^) + {p+p — > 
7r~)], the data are taken from (cf. and [|35| and references therein). The h~ from 



nucleon-nucleon interactions peak at midrapidity. The multiplicity increases with target 
mass, and the maximum of the distribution shifts to lower rapidities with increasing target 
mass. This is indicative of vr^ -pro duct ion after multiple collisions. For ?/ > 4, the yield of 
h~ produced in p-|-Au collisions appears to be lower than that for nucleon-nucleon and 



p-|-S collisions. Fig. 21 shows a comparison to minimum bias p+p, p-l-Ar and p-|-Xe data 



from experiment NA5 at the CERN-SPS [ p6|. 
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Figure 20. Rapidity distributions of negatively charged hadrons produced in minimum 
bias nucleon-nucleon, p-|-S and p-|-Au collisions. 
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Figure 21. Rapidity distributions of negatively charged hadrons produced in minimum 
bias p+p, p+S, p+Ar, p+Xe and p+Au collisions Note that the rapidity distribution 
for p+S collisions is higher than for p+Ar interactions due to different trigger conditions 
and event selections. 
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6.1.2. Nucleus-Gold Collisions 

The rapidity distributions of negatively charged hadrons produced in central d+Au, 
0+Au and S+Au collisions are shown in Fig. ^ The distributions exhibit a peak at y 
below Ucm = 3. The multiplicity increases with projectile mass, and the maximum of the 
distribution moves toward midrapidity with increasing projectile mass. 
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Figure 22. Rapidity distributions of negatively charged hadrons produced in central 
d+Au, 0+Au and S+Au collisions. 



6.1.3. Sulphur-Nucleus Collisions 

To investigate the target mass dependence for S-induced reactions, rapidity distribu- 
tions of negatively charged hadrons produced in central S+S, S+Ag and S+Au collisions 
are compared in Fig. |2^. The h~ distribution peaks at midrapidity for the symmetric 
S+S interactions. The multiplicity increases with the target mass, and the maximum of 
the distribution moves to lower rapidities with increasing target mass. Fig. ^ shows a 
comparison of the rapidity distributions of negatively charged hadrons produced in cen- 
tral S+S collisions and in minimum bias nucleon-nucleon interactions. The latter was 
scaled by the ratio of participant nucleons. The width of the distribution for S+S colli- 
sions is slightly narrower than the one for the nucleon-nucleon interactions. The mean 
negative hadron multiplicity per participant is higher in central S+S collisions than in 
N+N interactions by about 20%. 
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Figure 23. Rapidity distributions of negatively charged hadrons produced in central S- 
nucleus collisions. 
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Figure 24. Rapidity distributions of negatively charged hadrons produced in central S- 
S collisions and in minimum bias nucleon-nucleon interactions (scaled by the ratio of 
participant nucleons). 
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6.1.4. Extrapolation to Aw 

Since tlie rapidity distributions were measured over a wide kinematic range and cover 
almost all rapidities between target and beam rapidity, the spectra can be extrapolated to 
full phase space in order to obtain average total multiplicities. The results are compiled 
in Table |^. The extrapolation to target and beam rapidity was based on a gaussian 
fitted locally to the adjacent part of the rapidity distribution. The multiplicities increase 
with system mass. For comparison average total multiplicities reported from different 
experiments are also shown. For a review of multiplicities of negatively charged hadrons 
produced in nuclear collisions cf. [ |39| . 

The average multiplicity of negatively charged hadrons increases with the number of 
net protons. For not too asymmetric systems the ratio of the multiplicity of negatively 
charged hadrons per net (participant) baryon is shown in Fig. ^ This ratio is smaller in 
nucleon-nucleon and p+S interactions than in central S+S and S+Ag collisions. 
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Table 7 

Average multiplicity (in acceptance and total) of negatively charged hadrons produced 
in minimum bias nucleon-nucleon (N+N), minimum bias p-nucleus, deuteron-nucleus 
and central nucleus-nucleus collisions. The first column lists the reaction, the second 
column is the measured rapidity range and corresponding average multiplicity and the 
third column is the extrapolation to 4tt acceptance. 



< h > 


Reaction 


0.0 < y < 6.0 


4:71 


P+P [ < 


35 1 






2.85 + 0.03 


N+N [ 35] 


_ 


3.22 + 0.06 


p+Mg [ 


37 


1 




4.9 + 0.4 


p+S [ 


E 




5.0 + 0.2 


P+S 


5.7 ± 0.1 


5.9 + 0.2 


p+Ar [ 36 


1 


_ 


5.39 + 0.17 


P+Ag [ 


37 


1 




6.2 + 0.2 


p+Xe [ 


36 


1 




6.84 + 0.13 




8 < -u < 6 


An 


p+Au 


6 6 + 01 


9 9 + 9 

iy . iy _1_ 


p+Au [ 37| 




7.0 + 0.4 


p+Au 


38 






7.3 + 0.3 


'central' 
P+Au [ i 




9.6 + 0.2 




0.6 <y < 6.0 


Air 


d+Au 


17.2 ± 0.3 


23 + 4 




0.6 <y < 3.8 


An 


0+Au 


94 + 2 


137 + 9 


0+Au [ 


3J 






124.0 + 1.9 
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Air 


S+S 


90 + 1 


98 + 3 


S+S 


3 






94 + 5 


S+Ag 
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201 + 2 
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Figure 25. Average total multiplicity of negatively charged hadrons per net baryon pro- 
duced in minimum bias nucleon-nucleon, p+S, central S+S and S+Ag collisions as a 
function of the number of net baryons. 
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6.2. Transverse Momentum Spectra 

It is important to study the transverse degrees of freedom for produced hadrons in order 
to learn more about the degree of thermahzation in heavy ion colhsions and possible effects 
of collective flow on the spectra. Thermal model fits of the invariant cross sections (see 
Eqs. 1 and 2) as a function of Pt provide information on the production mechanisms for 
hadrons. It is known that pions, the majority of the negatively charged hadrons, are 
to a large part daughters from the decays of heavier resonances. In addition, possible 
transverse collective flow could further modify the spectra. It is thus not too surprising 
that the px spectra of negatively charged hadrons are, in general, incompatible with 
thermal parametrizations over the entire p^-range with the exception of p+S interactions 



[ 0, ^ Since the px distributions are found not to be thermal, a 

systematic analysis using Eq. 2 fits is excluded. Therefore, the data on transverse 
momentum distributions will be presented systematically without fits in the following 
sections. However, a fit to Eq. 2 was performed for all systems in the rapidity bin near 
midrapidity in a limited pr range (0.3 < pr < 1-0 GeV/c) in order to give a rough 
estimate of 'temperatures'. 

6.2.1. Proton— Nucleus collisions 

The transverse momentum spectra of negatively charged hadrons produced in minimum 
bias p+S and p+Au collisions are shown in Figs. |2^ a,b for four different rapidity intervals: 
0.8 < y < 2.0, 2.0 < y < 3.0, 3.0 < y < 4.0, A.O < y < 5.0. The corresponding spectra 
have been multiplied by the factors of 1000, 100, 10, and 1, respectively, for presentation 
on the same plot. 
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Figure 26. Transverse momentum distributions of negatively charged hadrons produced in 
p+S (a) and p+Au (b) interactions for the various rapidity intervals shown. The vertical 
scale is given in (GeV)"^. 
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6.2.2. Nucleus-Nucleus Collisions 

The transverse momentum spectra of negatively charged hadrons produced in central 
d+Au interactions and central 0+Au and S+Au collisions are shown in Figs. a-c for 
various rapidity intervals. The corresponding spectra were multiplied by factors of 1000, 
100, 10, and 1, respectively, in the case of d+Au collisions. For the heavy ion collisions the 
first two spectra (rapidities below midrapidity) were multiplied by factors of 100 and 10. 
The distributions especially at low rapidities (0.8 < y < 2.0) for central 0+Au and S+Au 
collisions clearly exhibit a deviation from a single exponential in pt with an enhancement 
at low transverse momenta in the backward rapidity region. 
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Figure 27. Transverse momentum distribution of negatively charged hadrons produced in 
central d+Au (a), central 0+Au (b) and central S+Au (c) collisions at different rapidities. 
The vertical scale is given in (GeV)^^. 
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The transverse momentum spectra of negatively charged hadrons produced in central 
S+S, S+Ag and S+Au collisions are shown in Figs. ^ a-c for various rapidity intervals. 
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Figure 28. Transverse momentum distribution of negatively charged hadrons produced in 
central S+S (a), S+Ag (b) and S+Au (c) collisions at different rapidities. The vertical 
scale is given in (GeV)~^. 
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6.3. Mean transverse momentum 

The invariant cross sections of negatively charged hadrons cannot be fitted over the en- 
tire pt range with a single exponential in or rriT- Thus, the mean transverse momentum 
is calculated for various systems and rapidity intervals, and the results are presented in 
Table |^. As an example, the rapidity dependence of (pr) in central S+A collisions is 
shown in Fig. ^ The mean transverse momentum reaches a maximum at midrapidity 
and exhibits little dependence on rapidity in the vicinity of y = 3. 

For y > 4.5, (pt) decreases rapidly. For the gold target (Fig. |3D|) there is a slight 
increase of {pt) in the target hemisphere with the projectile mass. For p+p interactions 
the (pt) is 340 ± 7 MeV/c for vr" [ |3|, ||] and 366 ± 2 MeV/c for all charged 

particles in 1.5 < y < 4.5 [ 
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Figure 29. Mean transverse momentum of negatively charged hadrons produced in central 
S+S, S+Ag and S+Au collisions as a function of rapidity. The open stars are S+S data 
reflected at y = 3. 
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Figure 30. Mean transverse momentum of negatively charged hadrons produced in d+Au, 
central 0+Au and S+Au collisions as a function of rapidity. 



Table 8 

Mean transverse momenta of negatively charged hadrons produced in p-nucleus, d- 
nucleus and central nucleus-nucleus collisions for various rapidity intervals. The sys- 
tematic error is less than 10 MeV/c. 



<Pt> [MeV/c] 


Reaction 


0.8 <y < 2.0 


2.0 <y< 3.0 


3.0 <y< 4.0 


p+S 


325 ± 6 


363 ± 8 


335 ± 8 


p+Au 


311 ± 5 


377 ± 6 


372 ± 7 


d+Au 


331 ± 4 


369 ± 5 


384 ± 5 




0.8 < y < 2.0 


2.0 <y< 3.0 


3.0 < y < 3.5 


0+Au 


353 ± 5 


373 ± 4 


374 ± 12 


s+s 


349 ± 9 


377 ± 4 


388 ± 11 


S+Ag 


346 ± 2 


376 ± 1 


385 ± 1 


S+Au 


359 ± 6 


384 ± 5 


371 ± 10 




4.5 < y < 5.0 


5.0 <y< 5.5 


5.5 < y < 6.0 


S+Au 


364 ± 2 


322 ± 3 


261 ± 4 
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6.4. Inverse slope parameters 

To facilitate comparison with other published results (see Table 10) fits to Eq. 2 were 
performed in the interval 0.3 < Pt < 1-0 GeV/c. The resulting inverse slope parameters 
T are listed in Table |^. A typical result of such a fit is shown in Fig. ^ which indicates an 
overall discrepancy between the concavely shaped data and the rather straight to convex 
curve corresponding to Eq. 2. The fit is thus carried out in the restricted range only, 
and this has always to be kept in mind in comparing to other results concerning h~ and 
7r~. Below 250 MeV/c the data points are all above the fit. This enhancement {low pt 



enhancement) is most pronounced at low rapidities in collisions with heavy targets [ ^ 



Fitting p+p data in a similar way and averaging the results from different experiments 



yields an inverse slope parameter at midrapidity of 142 ± 10 MeV [ pl^, 34, 40, HI 48 



Table 9 

Inverse slope parameters T2 (Eq. 2, fitted to the data in the interval 0.3 < pr < 1-0 
GeV/c) of negatively charged hadrons produced in p-nucleus, d-nucleus and central 
nucleus-nucleus collisions for various rapidity intervals. 



T2 [MeV] 


Reaction 


0.8<y< 2.0 


2.0 < y < 3.0 


3.0 <y < 4.0 


4.0 < y < 5.0 


p+S 


159 ± 13 


169 ± 13 


154 ± 10 


142 ± 12 


p+Au 


150 ± 8 


183 ± 10 


177 ± 9 


171 ± 10 


d+Au 


151 ± 9 


165 ± 10 


181 ± 12 


183 ± 11 




0.8 < y < 2.0 


2.0 <y< 3.0 


3.0 < y < 3.8 


4.2 < ?/ < 4.4 


0+Au 


189 ± 5 


188 ± 8 


175 ± 10 




S+S 


175 ± 4 


191 ± 7 


191 ± 6 


186 ± 5 


S+Ag 


184 ± 2 


199 ± 2 


201 ± 4 


185 ± 5 


S+Au 


182 ± 7 


200 ± 10 


190 ± 10 


187 ± 5 
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Table 10 

Inverse slope parameters T2 (Eq. 2, fitted to the data in the interval 0.3 < Pt < 1-0 
GeV/c) of negatively charged hadrons in comparison with other experiments. 



T2 [MeV] 



Reaction 



this paper 
h- 

0.8 <y < 2.0 



NA34 [ |3 
h- 

1.0 <y< 1.9 



p+Au/W 



150 ± 8 



161 ± 5 



0+Au/W 



189 ± 5 



167 ± 5 



S+Au/W 



182 ± 7 



165 ± 5 



NA44 [13, |3T| 



Reaction 



this paper 
h- 

3.0 <y < 3.8 



0.25 < mr - mo < 1 GeV/c^ 
3.0 <y < 4.1 



S+S 



191 ± 6 



148 ± 4 



S+Au/Pb 



190 ± 10 



156 ± 8 



S+Ag 




0,5 1 1,5 : 

Pt [ GeV/c ] 



Figure 31. Transverse momentum distribution of negatively charged hadrons produced 
in central S+Ag collisions. The solid line is a fit of Eq. 2 to the data in the interval 
0.3 < Pt < 1.0 GeV/c. The vertical scale is given in (GeV)~^. 
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The average transverse momentum of negatively charged hadrons exhibits only a slight 
increase with the system size, measured by the average total multiplicity of the collision, 
as shown in Fig. |32|, whereas the temperature of the net protons increases from about 150 
MeV in p-nucleus to 200 MeV in central 0+Au collisions and about 240 MeV in central 
S-nucleus collisions (Fig. ^). 



^ 400 




1 10 10 

<h'> 

Figure 32. Average transverse momentum in 2 < ?/ < 3 of negatively charged hadrons 
produced in p-nucleus, deuteron-nucleus and nucleus-nucleus collisions as a function of 
the average multiplicity. 
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^ 300 




Figure 33. 'Temperature' of net protons for p-nucleus, deuteron-nucleus and central 
nucleus-nucleus collisions as a function of the average multiplicity of negatively charged 
hadrons. 



7. Summary 

Wc have presented results of a systematic study, using the NA35 experiment at CERN, 
of rapidity and transverse momentum distributions of net protons and negatively charged 
hadrons. A large part of this data was analyzed employing a new automatic analysis 
technique developed by NA35. The most significant findings are: 

1) The shape of the rapidity distributions of net projectile protons {p — p) is similar for 
d+Au, central 0+Au and central S+Au collisions. Prom this we conclude that the stop- 
ping of projectile nucleons is independent of the size of the (light) projectile nucleus when 
incident on a heavy target nucleus (Au in this study). This suggests that there is no 
significant difference in the stopping mechanism as inferred from the proton rapidity for 
these light incident nuclei. 

2) The average rapidity loss {6y) for projectile protons which participate in the interaction 
is {6y) = 1.2 - 1.6 for various projectiles from protons to sulphur incident on a sulphur 
target and {6y) = 2.0 - 2.3 for a gold target. The {6y) are found to be approximately 
independent of the projectile mass for projectile masses ranging from protons to sulphur. 
Thus, the rapidity loss of projectile participant protons increases with the nuclear thick- 
ness of the target and is approximately independent of the nuclear thickness of the lighter 
projectile as might be expected. 

3) The number of net protons at midrapidity increases with projectile mass and with 
target mass. Since the mean rapidity shift measured for interacting projectile protons 
increases with target mass up to (Sy) = 2.0 - 2.3 for a gold target and since the rapidity 
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region for 200 GeV/nucleon incident energy is a total of yproj — Vtarg — 6, it is expected 
that the total number of protons at midrapidity (y = 3) will increase with both projectile 
and target mass. This results in a significant pile-up of protons at midrapidity for the 
heavier systems. 

4) The transverse momentum distributions for net protons {p — p) can be described by 
thermal distributions with 'temperatures' ranging from T 140 MeV for p + S near tar- 
get rapidity to T f5i 245 MeV for central S -|- S and S -|- Au collisions at midrapidity, using 
Eq. 1. Such inverse slope parameters are much higher than the Hagedorn temperature 
(T 150 MeV) and cannot be considered to be true thermal freeze-out temperatures. 
Our results may indicate the presence of radially ordered transverse flow in central S-l-A 
collisions. 

5) The average negative hadron multiplicity increases with both the target and projectile 
masses. The mean negative hadron multiplicity per participant is higher in central S-l-A 
collisions than in N+N interactions by about 20%. 

6) The mean transverse momentum of negatively charged hadrons at midrapidity is ap- 
proximately independent of the mass of the incident nuclei. There is a decrease in the 
mean transverse momentum when going away from midrapidity. This is observed for all 
colliding systems of this study. 

7) The spectra of negatively charged hadrons do not have a simple thermal shape. They 
exhibit deviations (above a thermal spectrum) at low {px < 250 McV/c) and high (px > 
1 GeV/c) transverse momenta. These significant deviations may be due to transverse 
collective flow and resonance production and decay as well as to the 'Cronin effect'. 

Acknowledgement 

We are very grateful to our scanning and measuring crews at Frankfurt, Munich 
and Warsaw for their careful and devoted work. This work was supported by the 
Bundesministerium fiir Forschung und Technologie and by the Leibniz Grant of the 
Deutsche Forschungsgemeinschaft, Germany, by the US Department of Energy (DE- 
AC03-76SF00098), by the Polish State Committee for Scientiflc Research (204369101), 
by the Polish-German Foundation (565/93) and by the Commission of European Com- 
munities (CI-0250YU(A)). 



46 



REFERENCES 

1. Proc. 11th Int. Conf. on Ultra-Relativistic Nucleus-Nucleus Collisions (Quark Matter 

1995) , Monterey, California, USA (1995), ed. by A.M. Poskanzer, J.W. Harris and 
L.S. Schroeder, Nucl. Phys. A590 (1995). 

2. Proc. 12th Int. Conf. on Ultra-Rclativistic Nucleus-Nucleus Collisions (Quark Matter 

1996) . Heidelberg, Germany (1996), ed. by P. Braun-Munzinger, H.J. Specht and R. 
Stock, Nucl. Phys. A610 (1997). 

3. J. Bachler et al., NA35 Collaboration, Phys. Rev. Lett. 72 (1994) 1419. 

4. H. Sorge, Z. Phys. C67 (1995) 479, C71 (1996) 178 (erratum). 

5. U. Ornik et al., Phys. Rev. C54 (1996) 1381; J. Bolz, U. Ornik and R.M. Weiner, 
Phys. Rev. C46 (1992) 2047. 

6. J. SoUfrank et al., Z. Phys. C61 (1994) 659; E. Schnedermann, J. SoUfrank and U. 
Heinz, Phys. Rev. C48 (1993) 2462. 

7. J. Bachler et al., NA35 Collaboration, Z. Phys. C52 (1991) 239. 

8. A. Bamberger et al., NA35 Collaboration, Z. Phys. C43 (1989) 25. 

9. J. Bartke et al., NA35 Collaboration, Z. Phys. C48 (1990) 191. 

10. J. Bachler et al., NA35 Collaboration, Z. Phys. C58 (1993) 367. 

11. T. Alber et al., Z. Phys. C64 (1994) 195. 

12. K. Werner, Phys. Rev. Lett. 62 (1989) 2460. 

13. J. Cleymans and H. Satz, Z. Phys. C57 (1993) 135. 

14. P. Braun-Munzinger et al., Phys. Lett. B365 (1996) 1. 

15. A. Leonidov, M. Nardi and H. Satz, Z. Phys. C74 (1997) 535. 

16. D. Brinkmann et al, Nucl. Inst, and Meth. A354 (1995) 419. 

17. S. Wenig, GSI-90-23 Report (1990). 

18. J.W. Harris et al., Nucl. Inst, and Meth. A315 (1992) 33. 

19. T. Alber et al., NA35 Collaboration, Phys. Lett. B366 (1996) 56. 

20. G. Roland, Univ. of Frankfurt PhD thesis (1992). 

21. B. Andersson, G. Gustafson and B. Nilsson-Almqvist, Nucl Phys. B281 (1987) 289. 

22. R. Brun et al., GEANT3 User Guide CERN/DD/EE/84-1 (1986). 

23. H. Appelshauser, Diploma thesis, Univ. of Frankfurt, GSI-93-09 Report (1993). 

24. L.P. Teitelbaum, PhD thesis, LBL-32812, Berkeley (1992). 

25. S.R. Tonse et al., NA35 Collaboration, Nucl Phys. A525 (1991) 689c. 

26. D. Rohrich, "Hadron Production in Nucleus-Nucleus Collisions at SPS- and AGS- 
Energies", pp. 505-516 in Hot and Dense Nuclear Matter, ed. by W. Greiner, H. 
Stocker and A. Gallmann, (NATO ASI Series B: Physics Vol. 335), New York and 
London: Plenum Press, 1994. 

27. J. Dodd, NA44 Collaboration, Nucl Phys. A590 (1995) 523c. 

28. F. Videbaek and O. Hansen, Phys. Rev. C52 (1995) 2684. 

29. A. Wroblewski, Acta Phys. Polonica B16 (1985) 379. 

30. A. Bialas, M. Bleszyhski, W. Czyz, Nucl Phys. Bill (1976) 461. 

31. l.G. Bearden et al., NA44 Collaboration, Phys. Rev. Lett. 78 (1997) 2080. 

32. A.M. Rossi et al., Nucl Phys. B84 (1975) 269. 



47 



33. T. Kafka et al., Phys. Rev. D16 (1977) 1261. 

34. M. Adamus et al., Z. Phys. C39 (1988) 311. 

35. M. Gazdzicki and O. Hansen, Nucl. Phys. A528 (1991) 754. 

36. C. De Marzo et al., Phys. Rev. D26 (1982) 1019. 

37. D.H. Brick et al., Phys. Rev. D39 (1989) 2484. 

38. J. Bachler et al., NA35 Collaboration, Z. Phys. C51 (1991) 157. 

39. M. Gazdzicki and D. Rohrich, Z. Phys. C65 (1995) 215. 

40. B. Alper et al., Nucl. Phys. BlOO (1975) 237. 

41. K. Guettler et al., Nucl. Phys. B116 (1976) 77. 

42. J.W. Cronin et al., Phys. Rev. Dll (1975) 3105. 

43. T. Akesson et al., HELIOS Collaboration, Z. Phys. C46 (1990) 361. 

44. H. B0ggild et al., NA44 Collaboration, Z. Phys. C69 (1996) 621. 

45. C. De Marzo et al., Phys. Rev. D29 (1984) 363. 

46. J. Whitmore, Phys. Rep. CIO (1974) 273. 

47. D. Rohrich et al., NA35 Collaboration, Nucl. Phys. A566 (1994) 35c. 

48. Y. Cho, Phys. Rev. Lett. 31 (1973) 413. 



